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Abstract 

In this paper we address the vector problem of a 2D half-plane interfacial crack loaded by a general 
asymmetric distribution of forces acting on its faces. It is shown that the general integral formula 
for the evaluation of stress intensity factors, as well as high-order terms, requires both symmetric 
and skew-symmetric weight function matrices. The symmetric weight function matrix is obtained 
via the solution of a Wiener-Hopf functional equation, whereas the derivation of the skew-symmetric 
weight function matrix requires the construction of the corresponding full field singular solution. 

The weight function matrices are then used in the perturbation analysis of a crack advancing 
quasi-statically along the interface between two dissimilar media. A general and rigorous asymptotic 
procedure is developed to compute the perturbations of stress intensity factors as well as high-order 
terms. 
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1 Introduction 



The theory of weight functions is fundamental in the evaluation of stress intensity factors for asymp- 
totic representations near non-regular boundaries such as cusps, wedges, cracks. The classical work of 
Bueckner defined weight functions for several types of cracks, both in 2D and 3D, as the stress intensity 
factors corresponding to the point force loads applied to the faces of the cracks placed in a homogeneous 
continuum (Bueckner, 1987, 1989). In particular, for a penny shaped crack in three dimensions, as well 
as a half-plane crack, Bueckner has introduced the weight functions corresponding to a general distri- 
bution of forces on the opposite faces of the crack, which also included the case of forces acting in the 
same direction. 

In the present paper, we use the term "symmetric" load for forces acting on opposite crack faces 
and in opposite directions, whereas the term "skew-symmetric" or "anti-symmetric" load will be used for 
forces acting on opposite crack faces but in the same direction, see Fig.[TJ Consequently, symmetric loads 
are always self-balanced, whereas skew-symmetric loads are not, in general (however, we shall consider 
only the case of self-balanced loading). Note that this terminology is different from that usually in use, 
where the words "symmetric" and "anti-symmetric" are intended in terms of the symmetry about the 
plane containing the crack, so that, for example, shear forces acting on opposite crack faces but in the 
same direction are considered to be symmetric (see for example Meade and Keer, 1984). 



crack 



crack 



symmetric component skew-symmetric component 

Fi gure 1: Symmetric and skew-symmetric parts of the loading in the case of single point forces. In practical cases, the 
loading is given as a combination of point forces and/or distributed forces in such a way it is self-balanced (see for example 
the problem analysed in Section [gj. 



Although in two dimensions skew-symmetric loading does not contribute to the stress intensity fac- 
tors, it becomes essential in three dimensions. Both symmetric and skew-symmetric loads on the faces 
of a half-plane crack in a three dimensional homogeneous elastic space are also considered by Meade and 
Keer (1984). 

The situation when the crack is placed at an interface between two dissimilar elastic media is sub- 
stantially different from the analogue corresponding to a homogeneous elastic space with a crack. Here, 
even in the case of two dimensions (plane strain or plane stress) , the stress components oscillate near the 
crack tip and also the skew-symmetric loads generate non-zero stress intensity factors. Also for the Mode 
III, where the stress components do not oscillate, there is a non-vanishing skew-symmetric component 
of the weight function. 

Symmetric twodimensional weight functions (i.e. stress intensity factors for symmetric opening loads 
applied on the crack faces) for interfacial cracks were analysed by Hutchinson, Mear and Rice (1987). 
The problem of symmetric three dimensional Bueckner's weight functions for interfacial cracks has been 
addressed by Lazarus and Leblond (1998) and Piccolroaz et al. (2007). On the basis of these results, 
Pindra et al. (2008) studied the evolution in time of the deformation of the front of a semi-infinite 3D 
interface crack propagating quasistatically in an infinite heterogeneous elastic body. 
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However, to our best knowledge, the skew-symmetric weight function for interfacial cracks has never 
been constructed for the plane strain case nor even for the Mode III deformation. On the other hand, 
in most appHcations, especially for the interfacial crack, the loading is not symmetric and consequently 
the effects of the skew-symmetric loading in the analysis of fracture propagation require a thoroughly 
investigation. 

The aim of this paper is to construct the aforementioned skew-symmetric weight functions for the 
twodimensional interfacial crack problem. To this purpose, we develop a general approach, which allows 
us to obtain the general weight functions, as used by WiUis and Movchan (1995), defined as non-trivial 
singular solutions of a boundary value problem for interfacial cracks with zero tractions on the crack 
faces but unbounded elastic energy. By taking the trace of these general weight functions on the real 
axis, one can arrive to the notion of Bueckner's weight function associated with the point force load on 
the crack faces (Piccolroaz et al., 2007). It is then shown that the skew-symmetric part of the loading 
contributes to the stress singularity at the crack edge and thus to the resulting stress intensity factors. 
These results are presented in Section [3] for the plane strain problem and in Section [7] for the Mode III 
problem. 

We summarize that the symmetric weight function matrix [[C/J for a plane strain interfacial crack 
has the form (compare with (|25)) in the main text): 

IUI{x^) = —^=\^B+^bA forxi>0, 

2rfoV27ra;i [ c{ J (1) 

[[{7Kxi)=0 forxiO, 

where i is the imaginary unit, B = , the superscript ^ denote transposition, and the bimaterial 

parameters e and do are defined in Appendix I A. II The quantities c^, cj~ appear as coefficients in the 
representation of oscillatory terms in the asymptotics of stress near the tip of the interfacial crack, as 
defined in the sequel (see ([30]) -(021)). 

The result for the skew-symmetric weight function matrix (U) , unpublished in the earlier literature, 
reads as follows (compare with (|26|) in the main text): 



{U){xi)^^m{^i) fora.i>0, 
/TT\( \ . oK-T*) /(-^i)21k L^kA f 



(2) 



where the bimaterial constant 7* and the Dundurs parameters a and c?, are defined in Appendix I A. II 
The conventional notations |/](a:^) = /(a;, 0+) — /(a;, 0~) and {f){x) ~ i(/(.T, 0+) + /(x, 0^)) are in use 
here to denote the symmetric and skew-symmetric components, respectively. 

Note that the result for the symmetrical weight functions is consistent with the representation for 
the stress intensity factors derived by Hutchinson, Mear and Rice (1987). 

For the Mode III case, the scalar symmetric and skew-symmetric weight functions are (compare with 
l(50l) and in the main text): 

-1/2 t 

tor xi > 0, ^ 77 



lUsj{x,)^{ V2^ ^ ' ', m = jlU,l (3) 

0, for xi < 0, 

where 77 is another bimaterial constant defined in Appendix lA.il 

In Section [4] we use both symmetric and skew-symmetric weight function matrices together with the 
Betti identity in order to derive the stress intensity factors for an interfacial crack subjected to a general 
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load. We summarize here that the integral formula for the computation of the complex stress intensity 
factor in the plane strain problem, K ~ Ki + iKu, has the form (compare with l|38p in the main text): 

K=-iMrMim / {m'^{x[-xi)R{p)ixi) + {U)^ix[-x,)RM{xi)}dxi, (4) 

where K ~ [K,K*]^, the superscript * denotes conjugation, jVti is defined in l(35|) and (|32| . R is 
defined in (|27|) and {p){xi), stand for the symmetric and skew-symmetric parts of the loading, 

respectively. 

For the computation of the Mode III stress intensity factor. Km, the integral formula has the form 
(compare with l(52|) in the main text): 

^iii = -/f / {(P3)(xi) + IIpbUxi)} {-x^r^'^dx^. (5) 

Needless to say, by replacing the loading with the Dirac delta function, these integral formulae 
immediately give the Bueckner's weight functions. Also, this approach allows us to derive the constants 
in the high-order terms, which are essential in the perturbation analysis (Lazarus and Leblond, 1998; 
Piccolroaz et al., 2007). 

Section[5]is devoted to the perturbation model of a crack advancing quasi-statically along the interface 
between two dissimilar media. We develop and prove an alternative approach consistent with the idea 
originally presented in WiUis and Movchan (1995). Both symmetric and skew-symmetric weight function 
matrices are essential in this perturbation analysis. 

A simple illustrative example is given in Section \E[ which shows quantitatively the influence of the 
skew-symmetric loading for interfacial cracks. 

Section [2] is introductory and contains auxiliary results concerning the physical flelds around the 
interfacial crack tip, which are necessary for the analysis presented in Sections [SHZl 

Finally, the Appendix lA.ll collects all bimaterial parameters involved in the analysis of the interfa- 
cial crack, including the well-known parameters as well as the new parameters appearing in the skew- 
symmetric problem. The parameters have been classifled according to whether they are involved in the 
symmetrical or skew-symmetrical weight functions. The Appendix IA.2I provides the theorems required 
for proving the approach presented in Section [H 



2 Preliminary results for an interfacial crack 

To introduce the main notations and the mathematical framework for our model, we start with the 
analysis of the displacement and stress fields around a semi-infinite plane crack located on an interface 
between two dissimilar isotropic elastic media, with elastic constants denoted by i^±,fi±, see Fig. [21 

The loading is given by tractions acting upon the crack faces. In terms of a Cartesian coordinate 
system with the origin at the crack tip, traction components are defined as follows 

afjixi,0^)^pf{xi), fovxi<0, J -1,2, (6) 

where p^{xi) are prescribed functions. 

The load is assumed to be self-balanced, so that its principal force and moment vectors are equal to 
zero. We also assume that the forces are applied outside a neighbourhood of the crack tip and vanish at 
infinity (specific assumptions on the behaviour of the loading will be discussed in the sequel). The body 
forces are assumed to be zero. It is convenient to use the notion of the symmetric and skew-symmetric 
parts of the loading 
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Figure 2: Geometry of the model. 



The solution to the physical problem is sought in the class of functions which vanish at infinity and 
have finite elastic energy. 

The results presented in Sections 12. II and [ 2.21 are known in the literature (see for example: Williams, 
1959; Erdogan, 1963; Rice, 1965; Willis, 1971), at least for symmetric loading. These results, together 
with the results on higher order terms presented in Section [231 which are less common in the literature, 
provide the formulae needed for the perturbation analysis developed in the sequel of the paper. 



2.1 Representation in terms of Mellin transforms 

In polar coordinates the Mellin transforms for the displacement vector and the stress tensor with locally 
bounded elastic energy are defined as follows 

/•OO /"OO 

u{s,6)^ u{r,ey-^dr, a{s,6) ^ (T{r,eydr, 







and they are represented by analytic functions in the strips — i?o < Re(s) < -doo and —70 < Re(s) < 700, 
respectively, where 'do,'&oo > (z?o + i?oo > 0), 7o,7oo > are constants related to the behaviour of the 
solution at the crack tip and at infinity, namely 

\ ■• I y Q(^j. »oo)^ r ~> cx), ^ ' [ 0[r T^=° r ^ 00, ^ ^ 

Correspondingly, the inverse transforms are 

uir,0)^— u{s,0)r-'ds, a{r,9) = — a{s,e)r-'-'ds, (9) 

Juii—ioo ^""^ Ju2—ioo 

where —i)q < oji < 1)00 and —70 < ^2 < 7oo- 



2.1.1 General solution to the field equations of linear elasticity 

In polar coordinates, with the origin at the crack tip, the Mellin transforms of stresses and displacements 
satisfy the relations (the material constants involved in these representations are reported in Appendix 

EH) 

CT± = Ct cos[(s + 1)61] + Ct cos[(.s - 1)61] + sm[{s + 1)6] + Cf sin[(s - 1)6*], 



6 



^ rr 



s - 1 



(1 - Zy±)CT(, ], Ug 



- 1 de 

\-v± d ^± 



where ct^ = at-^ 



Hence, we deduce 



cos[(s + 1)61] - cos[(s - 1)61] - ^-^^3=^ sin[(s + 1)61] - Cf sin[(s - 1)6*], 



s 

s_+ 1 

s 

1 



s + l 



— Cf sin[(s + 1)6*] - Cf sin[(s - 1)61] + cos[(s + l)6i] + Cf cos[(s - l)6i] 



[Cf cos[(s + 1)6*] + Cf cos[(s - 1)6*] + Cf sin[(s + 1)0] + Cf sin[(s - 1)61] 

+ ^^7TT^ [C± cos[(s + 1)0] + cf sm[(s + 1)^ 



C3±^cos[(,s + l)0] 



J cos[(s - 1)0] + [Ct sin[(.s + 1)0] - Cf cos[(s + 1)0]] 

The coefficients C^ depend on the loading and will be given in the next section. 
2.1.2 The boundary conditions and full field solution 

The coefficients C^ are obtained from the boundary conditions on the crack faces, namely 

c'eel^i^) ^P^ir), cr^eC^-.i^r) =q^{r), 

where p'^{r),q^{r) are prescribed functions given by (r) =pf[—r), q^{r) ~pf{~r), see Mishuris 
and Kuhn (2001) give the result in a compact form as the solution to the following system of algebraic 
equations 

Ctis) 



2 sinTTs 

s-1 



Cfis) J 



— sin TTS ± cos TTS 
zbCOSTTS sin TTS 



Cf{.s)+Cf{s) 



s + l 



— Cf(,s) + Cf(s) 



Dis)± 



D{s), 



(10) 



where D{s) = ^(s)F(s) and 
*(s) = 



cos TTS d* sin TTS 

— d* sin TTS cos TTS 



, F{s) 



Note that the prescribed boundary conditions appear in (fTO|) and in the definition of the vector F(s), 
in terms of the symmetric and skew-symmetric parts of the loading. The material constants d* and a 
are the Dundurs parameters, as described in the Appendix lA.il 
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Denoting the determinant of by S{s) — cos^ ns + d1 sin^ tts, we obtain the full field solution in 
terms of the Mellin transform as follows 



s — 1 1 
Ctis) = Trf*)cos7rs] + -IpI[{It d^)a cos ns] 



+ {q) [-(1 =F d*)d* sin tts] + -| 



/ 7 NT- /-, \ COS TTS 

(a, — ajct* smTTs + (1 =F a) 

SmTTS 



s-1 
2S{s) 



(p) [—(1 =F d,)c?, sin tts] + 



(d* — a)d^ sin TTS + (1 =F «)- 



sm TTS 



+ T d*) cos tts] + T d*)acos7rs] >, 



- — — <^ (p) [(1 + s ± (1 - s)d,) cos tts] + - H [(1 + s ± (1 - s)d, )a cos tts] 
Zd(s) 2 



+ {q) [(±(1 + s) + (1 - s)(i,)d, sin tts] 



COS TTS 

— (q(1 + s) + — s))d^ sin TTS — (1 =F — s) — ; 



1 



25{s) 



COS TT 5 

(a(l - s) + d^{l + s))d^ sin TTS + (1 T ")(! + s)— 



+ {q)[{l - s ± (1 + s)d*)cos7rs] + - s ± (1 + s)d^)a cos TTs]j. 

We note that some of the poles of the functions Cf{s) are obtained from the solution of the equation 

72 o;v,2 



S{s) ~ cos ITS + dl sin tts = 0, so that the poles are given by 



l-2n 



± ie. 



(11) 



where n is an integer, and e is the bimaterial constant defined in Appendix lA.il 

For the purpose of evaluation of residues, it will be useful to have at hand the derivative of S{s). 
This is given by S'{s) = 27r(— 1 + dl) sin tts cos tts, and its value at s = is as follows 



S' s 



±2TTid^ 



(12) 



This formula will be used in the next section for the derivation of the asymptotic estimates of stress and 
displacement fields near the crack tip. 



2.2 Asymptotic representations near the crack tip and the stress intensity 
factors 

The solution outlined above leads to the asymptotic representations of stress and displacement near the 
crack tip. Analysis of the singular terms in the stress components yields the complex stress intensity 
factor for the interfacial crack, where the stress singularity is accompanied by the oscillatory behaviour 
of the physical fields. 
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2.2.1 Asymptotics of stress and displacement near the crack tip 

Taking into account the assumptions on the appHed forces, the inspection of the results for cr(.s, 9) shows 
that o-{s,9) is analytic in the strip —1/2 < Re(.s) < 1/2. Therefore, in ^ we have 7o = 7oo = 1/2, and 
choosing W2 = 0, the inverse transform is given by 



0+ioo 



' — ioo 

By means of the Cauchy's residue theorem, we can get the asymptotics of cr(r, ^) as r ^ as follows 

pLU-{'ioo 



(T{r,t 



5^ Res] 



cr s, 



-s-l 



„±1 



1 

27ri 



<t(s,( 



-s-l 



ds, 



where lu < —1/2. The first two terms are the leading term, of the order 0(r~^/^), and the last term is a 
higher order term, of the order 0{r^),P > —1/2. The notation Res[/(s), s = s*] stands for the residue 
of / at the pole s = s*. 



Res[/(s), 5 = 5*] = 



1 



d"'-'fis)is-s*y 



(m- 1)! 



ds" 



where m is the order of the pole. 

It suffices now to compute the residues at the poles s = sf 



-1/2 ±ie, 



Res [a{s, 0) 



r ^ ^ , s 



= -1/2 ±ie] = 



■ Res 



\ s=-l/2±je 



lim 

-l/2±i( 
1 



-1/2 ±ie 
s - (-1/2 ±ie) 



27rd* 



where in the last equality we used the formula l|T2|). 

It can be shown that /3 = 0, so that the leading term asymptotics of the stress field is 



CT{r,e)=^_U2{r,e) + 0{l) 



(13) 



where 



= -l/2+ie 



+ [^(s,0)<5(s)],_i/2_r-V2+-} 



An inspection of the results for u{s, 0) shows that it is analytic in the same strip —1/2 < Re(s) < 1/2, 
except for a simple pole in s = 0. Since we seek the solution vanishing at infinity, the strip of analyticity 
for the function u{s, 6*) is < Rc(s) < 1/2. Therefore, in ^ we have t?o = 0, i?oo — loo — 1/2. Choosing 
ui in this interval and applying the Cauchy's residue theorem, we can write the leading term asymptotics 
of u{r, 9) as r ^ as follows 



-1 /•oj + ioo 

M(r, 9) = Res [m(s, 61) r-^ s = O] + ^ Res [m(s, 9) r-^ s = ] + — / it(s, 9)r-'ds 
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where w < —1/2. The first term of the order 0(1) corresponds to the translation of the crack tip and 
reads 



u{O,d)^Vo{0) = Q{0)wo, Q{9) 



cos ti sm fc 
— sin 9 cos ( 



(14) 



where 

Woi = Ml (0,0) 
1 



[1 - 2iy± ± 2d,iv± ~ 1)]7T / {p){r)dr + (-1 ± a)(i/± - 1) / Mir){logr)dr 



2TTfl± 
W02 =M2(0,0) 

= ^^|-[l-2z^±±2d,(z.±-l)]7r^ {q){r)dr + {-l±a){iy±-l)l^ MW(logOrf^| , 

are the Cartesian components of the translation of the crack tip. The second and third terms are of the 
order 0(r^/^), and the last term is a higher order term, which is of the order 0{rf^),(3 > 1/2. 

It can be shown that /? = 1, so that the computation of the residues leads to the asymptotics of the 
displacement field as follows 

M(r, 6) = Vo{e) + V,/2ir, 6) + 0(r), (15) 

where 

Note that the remainders 0(1) in lfT3|) and 0(r) in ifTSj) are rough estimates and will be refined in the 
Section [231 However, their physical meaning is already evident: the first corresponds to the so-called 
T-stress and the second to a rigid body rotation superimposed to a uniform deformation near the crack 
tip. 



2.2.2 The complex stress intensity factor K = Ki + iKn 

From the full field solution outlined in Sections 12.1.11 and 12.1.21 we obtain 



(s,0) + za.,(s,0) = - [{p){s) + i{q){s) + |H(s) + *|ra(s)} 



cos TTs + id* sin tts 



6{s) 



Applying the formula lfT3|). the leading term asymptotics may be written as 



(Jee{r,{)) + iare{r,Q) 



K 



-l/2+ie 



o(i). 



where 



-cosh(^6) {{p){s)+^{q){s) + |Ip1(.) ra(.)} 

TT L II) s=-l/2- 



is the complex stress intensity factor. 

Correspondingly, by means of the formula lfT5|) , we obtain 



(16) 



10 



where we used the notation [[/]]('') = /(^i ~ /('^i Note that the zero-order term Vo{0) present in 
ifTS]) disappears in the last formula. 

The formula for the complex stress intensity factor is then 

K = _y|cosh(^e)^°° {(p)(r) +*(<z)(r) + ||p](r) + (O} ^"'/'""d^- 

As expected, this representation is consistent with Hutchinson, Mear and Rice (1987) who considered 
a crack loaded by point forces a distance a behind the tip: 

p~^{r) = p^{r) = —P5{r — a), {r) = (r) = —QS{r — a), 

where S{-) is the Dirac delta function, so that 

{p){r)^-PS{r~a), MM=0, {q){r) ^ ~Qd{r ~ a), M{r) = 0. 

For this loading we obtain 

K = ■y/^cosh(7re) ^ {PS{r - a) + iQS{r - a)} r-^l'^'^^dr = ■y^cosh(7re)(P + iQ)a-^l'^-^\ 

which fully agrees with Hutchinson, Mear and Rice (1987). 



2.3 High-order asymptotics 

The asymptotic procedure involving the evaluation of stress intensity factors for cracks with a slightly 
perturbed front requires the high-order asymptotics of the displacement and stress fields near the crack 
edge. The procedure of the previous section can be extended to the high-order terms, and hence the 
high-order asymptotics of the stress field is given by 



<T(r, e) = 2_i/2(r, e) + T{e) + 3i/2(r, d) -f S{d)r -f ^^i^{v, 9) + Oir^), 
where (the superscript ^ denotes transposition, xi = rcosO, X2 = rsinO, Q{9) is given in (HJ 

1 



(17) 



T{e) = Res [<t(s, 6), s = -I] = TQ{0) 
S{9) = Res [a-{s,9),s^ -2] = 



1 =F a 



Q{0) 



X2 



X2 



dr 



Mir) 



dr 







dr 



Q ' (0), 



dr 



and 



, s = -l/2-l±ie] 



2TTd. 



{- [<^(.,0)<5(s)],=_,/2_i+,,r'/2- + [ais,e)S{s)l_,/,^,_^y/'+^^}, / = 1,3. 



Note that T is the T-stress given by 



T = 



1 T a 



r) 



dr 
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Applying the formula (fT7|) to (|16p . we finally obtain 

^ee{r. 0) + *c7.,(r, 0) = J^r-'/^+^^ + -^r''^+^^ + ^r3/2+- + 0{r% 
v27r v27r v27r 



where the constants in the high-order terms are 

B = -^coM^e)j^ [{p){r)+i{q){r) + '^Ip^^r) + z^lql{r)] r-^'^-^^dr. 
If the loading is smooth enough, then we can integrate by parts to obtain 

" - -{^^WT^wkir)[ {<!'>"(^) +*<">"(••) + fM-c-) M"(-)} ^-'"-"'i- 

Correspondingly, we obtain for the high-order asymptotics of the displacement field the expression 
u{r, 9) = Vo W + Vi/2(r, 9) + Vi{9)r + V^/^ir, 9) + V2{9y + V5/2{r, 9) + 0(r3), 

where 

Vi{9) = Res [m(s, 61), s = -I] = Q{9)r-^wi{r, 9), 1^1,2, 

yi/2 =^Rcs[u{s,9)r-'\s ^ -l/2±ie\ 
± 



and 



wii = -^^" f 
27r/i± 

1 =F a 



Wl2 



27r/i± 



r 


/"OO 




r 


/•OC 

Jo 


r 



W2l^^^^l->^±)xrX2J^ Mir)^ + [i4+4)-{3-2,^±){xl-xl)]J^ MM^j, 

1 =C n- r r°° dr f°° dr^ 

W22 - ^ + + (1 - 2z.±)(.T? - xl)] M{r)-^ + Ay±x^X2 M[r)-^^ . 

We note that wi{r, 9) is homogeneous of degree I in the variable r and that wn and wi2 are the Cartesian 
components of the local rigid body rotation superimposed to a uniform deformation near the crack tip 
mentioned above. 
Hence 



l/2-l-ieV2^ 3/2 + ieV2^ 5/2 + ^6^2^ 
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3 Symmetric and skew-symmetric weight functions for interfa- 
cial cracks 



In this section, we introduce a special type of singular solutions of a homogeneous problem for an 
interfacial crack. The traces of these functions on the plane containing the crack are known as the 
weight functions, and they are used in the evaluation of the stress intensity factors in models of linear 
fracture mechanics. The notations [L''] and (U) will be used for the symmetric and skew-symmetric 
weight function matrices, as outlined in the sequel of the paper. 

We refer to the earher publications by Antipov (1999), Bercial et al. (2005) and Piccolroaz et al. 
(2007) for the detailed discussion of the theory of weight functions and related functional equations of 
the Wiener-Hopf type. Here we give an outHne, required for our purpose of evaluation of the weight 
function matrices for interfacial cracks. 

The special singular solutions for the interfacial crack are defined as the solutions of the elasticity 
problem where the crack is placed along the positive semi-axis, xi > 0, the boundary conditions are 
homogeneous (traction-free crack faces) and satisfying special homogeneity properties. In particular, 

(a) the singular solution U ~ [Ui, 1/2]^ satisfies the equation of equilibrium; 

(b) lUj = when xi < 0; 

(c) the associated traction vector acting on the plane containing the crack, S — [E2i,S22]^, is con- 
tinuous and S = when a;2 = and a;i > (homogeneous boundary conditions); 

(d) [/ is a linear combination of homogeneous functions of degree —l/2 + ie and —1/2 — ie. 

We enphasise here that the domain for weight functions (where the crack is placed along the positive 
semi-axis) is different from the domain for the physical solution (where the crack is placed along the 
negative semi-axis). The reason for this is to have the fundamental integral identity, introduced in the 
sequel, eq. ((28l) . in convolution form (see Piccolroaz et al., 2007, for details). 

3.1 The Wiener-Hopf equation 

3.1.1 Factorization and solution 

Let us introduce the Fourier transforms of the crack-opening singular displacements and the correspond- 
ing traction components 



The superscript indicates that the function |i[/]]"'"(/3) is analytic in the upper half-plane C+ = {/? £ 
C : Re/3 G (— cxd, cxd), Im/3 € (0, 00)}, whereas the superscript ~ indicates that the function S (/3) is 
analytic in the lower half-plane C~ = {/3 S C : Re/3 € (— 00, 00), Im/3 G (— oo,0)}. 

These functions are related via the functional equation of the Wiener-Hopf type (Antipov, 1999): 




4g(/3)S (/3), /3eM, 



(18) 



where 6 = (1 — + (1 — i^_)//i_ and 



1 



i sign(/3)c?* 



isign(/3)d* 



1 
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The factorization of — 6/|/3|G is given by 



(19) 



where 

X+{(])^do 



cosi?+ — sini?+ 
sini?+ cosB'^ 



cos B — sin B 
sin B^ cos i?^ 



do — {I — d*)^^**, and the limit values of the functions , as Ini(/3) — > 0''', are 



,1/2 



and 



/31/2, > 

i(-/3)i/2, /3<0 



-elog(Tz/3), /3eC±. 

From the definition of singular solutions given in the introduction of this section, it follows that the 
asymptotic behaviour of the weight functions is given by 

[t/Kxi)~x7'/'Fi(xi), XI ^0+, 

- (-Xi)-3/2F2(2;i), XI ^0", 

lUKxi) x^^^^Fsixi), xi^co, 
S(xi) - (-a:;i)-3/2F4(a:;i), -c», 



where Fi, F2, -F3, F4 are bounded functions. The appHcation of the Abelian type theorem I A. II and the 
Tauberian type theorem IA.2I (see Appendix IA.2P gives 





/3eC+, 


/3- 


+ 00 


S"(/3)^/3y'F2(/3), 






^ 00 


[[/] + (/?) ^/3;^/'F3(/3), 


(3 e C+, 


P- 




S(/3) ~ p]/'F4{f3), 






-^0, 



where Fi, F2, F3, F4 are bounded functions. Substituting the representation (fT9| for the matrix 
— 6/|/3|G(/3) into the Wiener-Hopf equation ifTSj) . we obtain 



(20) 



Since [X'''(/3)]^^ and [X~{j3)]~^ are bounded functions and taking into account the asymptotic be- 
haviour of the weight functions at infinity, it follows that the LHS and RHS of l(20|) behave as 0(1), 
/3 00, and from the Liouville theorem, they are equal to the same constant e = [ei, 62]^, where ei, 62 
are arbitrary, so that 



[[/]+(/?) = -ijX+(/3) 



ei 

62 



f3eC+, E (/3) = -i^x-(/3) 



ei 

62 



(21) 
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3.1.2 The basis in the space of singular solutions 



Prom now on, we shall use the notion of weight functions, defined as traces of the singular displacement 
fields obtained in the above section. It is evident from the solution obtained above that the space of 
singular solutions is a 2-dimensional Hnear space. Two Hnearly independent weight functions can be 
selected as follows. Setting ei = 1, 62 = in ((2T|) . we get the first weight function as: 



1/2 



cos = 



It^2r = TrksinB+ = :^/3+ 



Zj2i — 



^22 — 



1/2 

+ 

1/2 



dOr,-l/2 

2 



eoP4 



1 

eo' 
1 



-A 



+ 



eo ^ 



1/2 



■ cos B = — - 



1 



^3^/2 ( + e^pz' 



P. 



bdn 



■ ainB 



2bdo 



1 

eo 



-/?!! - eo/?: 



(22) 



where eo = e" 

Setting ei = 0, 62 = 1 in l(2T|) . we get the second weight function as: 



/?+/^ ^ V eo 

_ cos = y /?+ (^eo/?;^ + -/?+ 

-2- P'/ 

= . . sini3 = 



bdo 



2bdo 



1 

eo 



' /3y' (-/?!!- eo/3: 



^22 — ^ 



/3 



1/2 



1 



6do 



cosS- = -— — /3i/2 ( — /3!_^ + eoPZ 



2bd, 



1 



eo 



(23) 



In order to use a compact notation in the sequel of the paper, we collect the weight function components 
together with components of the corresponding tractions in matrices as follows 



s = 



^21 ^21 



■^22 



■^22 



Note that the Hmit values of the functions /3^"^, as Im(/3) 0+, are 

r r , /3 > 



/?<0, 



(i7 



eoi 



/3>0, 
/?<0, 



whereas the limit values of the functions /3_", as Im(/3) — > , are 



/3>0, 
{-PY'el /3<0, 



/3>0, 
/? < 0. 
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It is clear that the choice of a basis of two linearly independent singular solutions is not unique. 
Bercial-Velez et al. (2005) provided another set of linearly independent weight functions, namely 



pill 



pill 



1/2 



^0 



^1 
c7 



eoct 



/3 e 



It can be shown that these weight functions are linear combinations of (|22| and l(23|) . 



P2I 



z(c++ci) 



3.2 The half-plane problem and full representation of weight functions 

In this section, we will construct the full representation of singular solutions for the whole plane. This 
is needed in order to compute the skew-symmetric weight function matrix (U). The complete singular 
solutions can be constructed by solving a boundary value problem for a semi-infinite half-plane subjected 
to traction boundary conditions at its boundary. 

Let us consider first the lower half-plane. A plane strain elasticity problem is usually solved by means 
of the Airy function. Introducing the Fourier transform with respect to the xi variable, the problem is 
most easily solved considering the stress component CT22 as the primary unknown function, so that the 
problem reduces to the following ordinary differential equation 

- 2/324'2 + - 0, 

where a prime denotes derivative with respect to X2 ■ The general solution is then 



^22 



{P,X2) = {A2+X2B2)e^^^''^ an{P,X2)^--^a'^2, a2i{P,X2) = 



JjO-22, Cr2l(/3,X2) = -■^cr22> 



Correspondingly, the Fourier transforms of the displacement components are 



Ui{P,X2) 



-{0-22 - (1 - i^-)o-o}, U2iP,X2) 



1 



r{a^2 + (l-'^-K}, PeR, 



where ctq = o^n +^22. The boundary conditions along the boundary X2 = 0^ are defined by 

^22(/3,a:^2=0-) = E2"2(/3), ^2i{P, X2 ^ 0') ^^2i{f3), (3&R, 
where S^jj ^21 the tractions along the interface derived in the previous section. It follows that 



a22(Ax2=0-)=A2 = S22(/3), a2i(/3,X2=0-) = --(A2|/3|+B2) = S2i(/3), P€R, 



and thus 



A. 



B2 ^ ip^2i - \P\^22, 
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The full representations of Fourier transforms (with respect to xi) of the required singular displacements 
and the corresponding components of stress are 

a22(/3,a;2) = {tPx2^2i + (1 - |/?k2)E2"2}el'5|-^ 

aii(/3,X2) = {-i{2sign{f3)+pX2)^2i + (1 + \f3\x2)^22}e^''^^' , 

a2iif3,X2) = {(1 + 1/31x2 +i/3x2E2~2}el''l-^ 



Ui{P,X2) = 
U2{P,X2) = 



1 



2/i_ 
1 

2^ 



X2 



2(1 



sign(/3)x2 



^21 +i 
1 - 2i^_ 



sign(/3)a-2 + 



1 - 2z/_ 



■^21 



2(1 



X2 



^22 



.1/3^2 



= 1/3^2 



^_ with /i+ 



For the upper half-plane, we find the same equations, subject to replacing \(3\ with 
and ly- with 

3.3 New results for skew-symmetric weight functions 
3.3.1 The Fourier transform representations 

It is possible now to derive the jump and average of Fourier transforms of the singular displacement 
functions across the plane containing the crack. The traces on the plane containing the crack are given 

by 

i(l-2i.+ ) 



Ui{(3)=ui{P,X2 = 0+) 

U2{P)=U2{f3,X2=Q+) 

Ui{P)=ui{f3,X2 = Q-) 
C/2(/3) =W2(/3,X2 =0-) 



i{l - 2v+) \-v+ 
i{l - 2i/_) 1 - zy_ 



^21 



22 
^22 



so that, we obtain in matrix form 



i sign(/3)d* 



2/x_/3 ' /i_|/3| 



-i sign(/3)(i, 
1 



^21 

^22 
^21 

^22 



Zj2i 


^^21 


-^22 


-^22 



and 



(ul) {ul) 
(ul) (ul) 



ba 



-i sign(/3)7* 
1 



^^21 ^^21 



where 7* is a material parameter, 



i sign(/5)7* 
'~ 2/._(l-j.+ )-2/i+(l-j._)- 



^22 ^22 



/3eR, 



/3 e R, 
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Note that the above equations can be rewritten as 



b 


" 1 ■ 


s 




" -1 " 


W\ 


. 1 . 






1 



(U) 



ah 



1 
1 



-1 

1 



/3 e K, 



so that we can derive the decomposition of {U) in the sum of " + " and " — " functions as follows 

-1 



{U)^-lUr + a{d.-l*)Yp 



1 



(24) 



where the last term on the right-hand side is a " — " function. 
3.3.2 The weight functions - Fourier inversion 

After the inversion of the corresponding Fourier transforms, the weight functions, which will be needed 
for the computation of stress intensity factors, are as follows. The symmetric weight function matrix 
[[J7]](a;i) is equal to for xi < 0, whereas for xi > it is given by 



lUl}ix,) = 



Xl 


1/2 


fx^' 




2do 






Cl 


ix-^ 


1/2 


(xT' 


xf 


2do 




Kct 





(25) 



xi) = iun{x,). 

The skew-symmetric weight function matrix (U){xi) is equal to — [[t/K^i) for > 0, whereas for xi < 



it is given by 



(uDixi) = 

(C/])(.Tl) = 



m(d, -7.)(-xi)-i/2 



MlV2^ 



i-xiY 



i-xiY 



a{d^ - 7*)(-xi) 



4:dlV2n 



hx,Y 



(26) 



(C/f)(.Tl) = -(C/])(.Tl), 

mix,) ^ (uDix^). 

The function S(xi) is equal to for xi > 0, whereas for xi < it is given by 

(-Xi)-3/2 



T.l2{xi) 



2bdlV2^ 

z(-Xl)-3/2 

2bdlV2n 

^li(xi) = -S22(a;i), 

^Uxi) = ^liixi). 



1/2 + ie 



1/2 -ie 



i-XiY 
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It will be shown in the sequel of the text that we only need the inverse transform of [[[/]]"'" (see l(24|) ) 
in order to compute the stress intensity factors for both symmetric and skew-symmetric loading. 



4 The Betti identity and evaluation of the stress intensity factors 

Here we develop a general procedure for the evaluation of coefHcients in the asymptotics of the stress 
components near the crack tip. This includes the stress intensity factors as well as high-order asymptotics. 
In particular, the coefficients near the higher order terms require appropriate weight functions which are 
shown to be derived via difi^erentiation of the weight functions [[[/]] and {U) along the crack. 

4.1 The fundamental Betti identity and its equivalent representation in 
terms of Fourier transforms 

In this section, the notations u = [ui, ^2]^ and cr = [(T21, 0-22]^ will be used for the physical displacement 
and traction fields discussed in Section [51 The notations U = \Ui, U2\^ and S = [£21,222]^ will apply 
to the auxiliary singular displacements and the corresponding tractions obtained in Section 13.3.21 We 
note that U is discontinuous along the positive semi-axis xi > 0, whereas u is discontinuous along the 
negative semi-axis xi < 0. 

Similar to WiUis and Movchan (1995) and Piccolroaz et al. (2007), we can apply the Betti formula 
to the physical fields and to the weight functions in order to evaluate the coefficients in the asymptotics 
near the crack tip. In particular, applying the Betti identity to the upper half-plane and lower half-plane, 
we obtain 



/(X2=0+) 

and 



{[/^(x'l -a;i,0+)i?cr(a;i,0+) - Yi^ {x[ ~ xi,Q+)Ru{xiM)} dxi = 0, 



(2:2=0) 



/ \u^{x^ - xi,Q-)R(t{xi,0-) - 11^ (x\ - xi,0~)Ru{xi,Q^) \ dxi = 0, 

respectively, where 

?) 

is a rotation matrix. Then, subtracting one from the other, we arrive at 
|(7^(x'i -xi,0+)J?cr(xi,0+) - U'^ {x[- xi,Q-)RcT{xi,Q-) 

- S^(x'i -xi,0+)fiM(xi,0+) - S^(a;'i -xi,0")ilM(cci,0") ^ dxi = 0. 
The traction components acting on the xi axis can be written as 

(t(xi,0+) =p+(.Ti) + crW(xi,0), (t{xi,0~) = p-{xi) + a^+\xi,0), 

where p'^i^xi) = a(xi,0^)H{—xi), p^{xi) ~ <t{xi,0^)H{—xi) is the loading acting on the upper and 
lower crack faces, respectively, and cr(+) (xi,0) is the traction field ahead of the crack tip, with H{xi) 
being the unit-step Heaviside function. Consequently, the integral identity can be written as 

/ {U^{x'i - xi,{)+)Rp+ {xi) + U^{x[ ~ xi,0+)R(T^+\xi,0) 

J{X2=0) 

-U^ix[ - xi,0-)Rp-{xi) - U'^ix'i - xi,0-)R(t'^+\xi,0) 

- 'E'^{x[- xi,0+)Ru{xi,0+)-'E'^{x[- xi,0-)Ru{xi,0-) jdxi = 0. 
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From the continuity of cr(+) and S we get 

/ f [i71 ^ (x'l ~x,)Ra^+\x,,0)^-E'^{x[^x,,0)Rlul{xAdx,= 

- / {x[ - xi,0+)Rp+{xi) - {x[ - xi,0-)Rp- {xi)]dxi. 

Introducing the symmetric and skew-symmetric parts of the loading |[7]) , we finally obtain 
/ -xi)i2(T(+)(xi,0) - {x\- xx,^)R\u\{x{)\dx^ = 



'{X2=0} 

'(X2=0) 

Taking the Fourier transform of (|28|) we derive 
IU}+^ Ra+ Rluj- = 



(28) 



lUl+^R{p)-{U)'Rlpl f3e 



(29) 



4.2 Stress intensity factors and high-order terms coefficients 

In this section, we will derive the stress intensity factors as well as high-order terms from the Betti 
formula (|29| obtained in the previous section. The complex stress intensity factor will be denoted by 
K = Ki + iKii, and analogous notations will be used for high-order terms: A = Ai + iAu, B = Bi + iBu. 
The asymptotics of the physical traction field as xi 0+ are as follows 

-1/2 1/2 3/2 

X^ , X^' , . I-/ 5/2, 



^W(^-i) - ^S{x^)K + -2^5(xi)A + -^S{x^)B + 0(x?/^), 
I'J I'K iy 111 Zy III 



where 



S{xi) 



and K = [K,K*]^, A = [A, ^*]^, B ~ [B,B*]^ , where the superscript * denotes conjugation. The 
corresponding Fourier transforms, as /? ^ oo, /3 G C+, are (compare with Piccolroaz et al., 2007) 



^-1/2 ^-1/2 ^-1/2 



4/3 



4/32 



T3(/3)B + 0(r 



where 



(i + «)V^ 



(l-z)V^ 



i = 1,2,3, 



(i + i)V^ 

'2r(5/2±ie)' 



1 2r(l/2±ze)' 2 2r(3/2±ie)' ^ 

and r(-) denotes the Gamma function. 

The components of the displacement physical field as xi 0~ have the form 



(30) 



(31) 



(32) 



M(a;i) = ^{-x,f''-U,{x,)K ^{-x,f/^U2ix^)A - 
v27r \/2tt 



/27r 



M~xif/^U3ixi)B + 0[i~xi\ 



7/21 
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where 



l(-Xi)" 



2j -1 + 2ie 2j -1- 2ie 



J = 1,2,3. 



L 2j -1 + 2ie 2j -1- 2ie 
The corresponding Fourier transforms, as /3 cx), /3 e are (compare with Piccoh-oaz et al, 2007) 



where 













c-eo 






4 


c-eo . 



/3:^'^V3(/3)B + 0(/3-'/'), (33) 



J = 1,2,3. 



From these asymptotics we can derive estimates of the terms in the LHS of the Betti identity (|29l 

as /3 ^ oo. 



/3e C^ 
/?e C- 



where 



do 


-s 




^4^7 







j = 1,2,3. 



(34) 



(35) 



Let us introduce the following notation for the RHS of the Betti identity l|29p 

*(/?) - -pl+^(/3)H(p)(/?) - {U)^{(3)RmW. /? e R- 

Note that for any possible loading, ^{(3) — *■ 0, as /? ^ ±oo. In fact, if the loading is given by point 
forces in terms of the Dirac delta function, then '4'(/3) ~ 0(|/3|~^/^). If, for example, € Li{R) or 
(p±)' G Li(R) then *(/3) = o(|/3|-i/2) or *(/3) = o{\f3\~^/^), respectively. 

Consequently, we can split ^^{(3) in the sum of a plus function and a minus function 



where 



*±(/3) 



27ri 



Moreover 'i'^(/3) ~ /3 cx), /3 e C^, and, from the estimates we can conclude that the 
inhomogeneous Wiener-Hopf equation (|29|) has the unique solution 



lU}+^RW+ = ^+, /3eC+ and S i?[[Mr = e C". 

From these identities we can extract asymptotic estimates. For this reason, let us first assume that 
G ^(M), where S{R) is the Schwartz space of rapidly decreasing functions (Zemanian, 1968), so that 
*(/3) decays exponentially as /3 — > ±oo. Then we have, for /3 ^ oo, /3 G C^, 

*(t) 



1 /•= 



t~ (3 ' 2TTi 



^{t)dt-(3-^ / t^{t)dt + 0{(3~^). 



(36) 
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Comparing corresponding terms in l|36]) and (|34|) . we obtain the following formulae for the complex 
stress intensity factor and the coefficient in the second-order term 



respectively. 
Note that 



^ r°° 

K = ^>tr' J {iur^m{p){t) + {u)^{t)Rmit)} dt, 

A = ^^M^' J t{lUi+^{t)R{p){t) + {U)^{t)Rlpm}dt, 

1 /-oo 

= lim — / {lUr^{t)R{p){t) + {U)^{t)Rmit)}e-"^'dt 

= lim r {lU}'^{x[-x,)R{p){x,) + {U)^{x[-x,)RM{x^)}dxi. 



The formula for the complex stress intensity factor becomes 

rO 



K = -tM^'\im f {lUVix[^x,)R{p)ixi) + {U)^ix\-x,)RMixi)}dx, 

^l-*0 J-oo 

Similarly, we can develop the integral in l(37|) 9 to obtain 
1 /-oo 

■V r fun^f ' Myd{p)ixi) , ,T, , \DdlpKxi)\ 

= ilim / UUl (x^~xi)R — h (J7) [x^-xi)R — } dxi. 

<^oJ_^ I dxi dxi J 



The formula for the coefficient in the second order term is 



(37) 



(38) 



A^M^Uim f \im^(x', - x^)r '^^^]^''''^ + ([/)^(x; - x,)R '^^^f'''-^ \ dx,. (39) 
^'i^"J-oo I dxi dxi ) 

If the loading p^ is not smooth enough, singular or it is a distribution (such as the Dirac delta 
function), then we can replace p"^ with a sequence {p^} such that p^ € iS(M) and p^ p^ in 
<S(M), when n ^ oo. Here 5(M) denotes the space of distributions (Zemanian, 1968). Therefore, the 
representations ((37|) . ((38|) . ((39|) can be interpreted in the sense of distributions and we can integrate 
formally l(39|) by parts to obtain 

1,. (dlur(x\-xi)„, , d(Uy(x[~xi) , 

A = A<2"Mmi / : ^■R(p)(xi)+ ^ ' \ ^ ^-^RM{xi)\dxi. 

^'i^oj-oo I dxi dxi ) 



5 Perturbation of the crack front 

For the 2D case, the perturbation considered here is related to an advance of the crack front by a distance 
a. We denote quantities relative to the original unperturbed crack problem by subscript 0, and quantities 
relative to the perturbed crack problem by subscript *. 
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5.1 General settings 

The complex stress intensity factor for the original unperturbed problem is given by 

Ko = -zMrMim / {m'^{x[-xi)R{p)ixi) + {U)^{x\~xi)RM{xi)}dx,. 

For the perturbed problem, if the coordinate system is shifted by the quantity a to the new crack tip 
position, the complex stress intensity factor is given by the formula 

.0 

K^{a) = -iM^^ lim / {fU}'^ {x[ ^ xi)R{p){xi + a) + {U)^ {x[ ^ xi)RM{xi + a)} dxi. 

^'l^O J-oo 

The computation of the first order variation of the complex stress intensity factor is now straightfor- 
ward 



da 



a=0 

+a) - {p){xi) 



tMrMim lim / huf{x[~xi)R- 



-zMrMim r {lUVix[-x,)R^MM + (^u)\x[-x,)R^M(^\dx, 
x'i~*oJ^rx, I dxi dxi ) 



Similarly, for the second order term, we obtain 
dA^{a) 



da 



= AI9 ^ lim 



hm / <^ |(7J| [Xi - xi)R -7-^ h {U) [x^ - xi)R — ^ > dxi. 

x[~^o J t ax^ ax-^ J 



a=0 

We assume here that all the derivatives exist, otherwise, as previously, the formula should be under- 
stood in the generalized sense. 

5.2 Computation of the perturbations of stress intensity factor according to 
the approach of Willis and Movchan (1995) 

The Betti formula for the original crack problem and the perturbed crack problem reads 

lur^Rwi^s Riuu^-iur^R{p)-{u)^mpi /Se^, 

respectively. According to l(30|) and l(33|) . the asymptotics of physical fields in the unperturbed problem 
are 

.-1/2 .-1/2 

2 2 

For the physical fields in the perturbed problem, crl and ul, we can write the respective transforms 
in the new coordinate system related to the new position of the crack tip 



Wi{l3,a) = / ai{xi)e'''^'dxi = e'^'^ / (r,{y)e'^ydy = e'^'^Wt if3, a), 
lullif3,a) = r lMi]](xi)e'^-^dxi = e^^^ [ iu4iy)e^0ydy = e'^l«.r a)- 
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Note that and |it*]l~ are " + " and " — " functions, respectively, with the same asymptotics as in (|4T|) 
subject to replacing Kq with -Rr*(a) and Aq with A^(a). 
Subtracting (|40)) 9 from (|40|) i . we obtain 



= 0, /3g 



(42) 



Note that for any a > 0, the function e^'^°-wj^ is a " + " function which decays exponentially as /? ^ oo, 
/3 £ C"*", whereas e'''°[Mj,]~ is a " — " function which grows exponentially as /3 oo, /3 € C~. 
Following WilHs and Movchan (1995), we expand the exponential term as a ^ 0+, so that 



S Riluoj 



0, 13 e 



The substitution of two-terms asymptotics for the unperturbed physical fields, ctq, fuo^ , and for 
the perturbed physical fields, leads to 



-{Mi{Ko - K^{a))l3Z^ - iaMiK^{a) - la, 



M2A,{a)pz^} + 0{p-'') = 



{Mi{K„ - K,{a)) - iaM2A,{a)}{(3+^ - pZ^) + 0{f3-^) 



We get 



AK = K^{a) - - -mMr^AisAo = a 



1/2 + ie 
1/2 -it _ 



(43) 



and finally 



AX - a(l/2 + ie)^o, AK* ^ a{l/2 ~ ie)Al, a ^ 0^ 



Note that it is assumed here that -FC*(a) and Af,{a) are continuous in a = 0. This fact will be proven 
in the next section. 



5.3 An alternative computation of the perturbations of stress intensity fac- 
tors and high-order terms 

Here we provide a rigorous asymptotic procedure leading to the computation of perturbations in the 
stress intensity factors as well as high-order terms. The Abelian and Tauberian type theorems outlined 
in Appendix IA.2I will be used for this purpose. 
We rewrite now the equation l(42|) as follows 

(44) 



The function e~^^°-T, iljito]]^ is a " — " function which decays exponentially as /3 ^ oo, /3 e C^, 
whereas e~^^"'\U\'^"^ Ra^ is not a proper " + " function, because of the exponential growth as /3 ^ oo, 
/3 G C"*". However, for /? G K, e~^^°-\U\^"^ Ra^ ~ 1//3, as /3 — > ±cx), so that it can be decomposed as 



e-''P<'lU^+^ Ra+ = a+ -a- 



where 



1 

27ri 



dt, (3 e 
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Now, from the solution of the Wiener-Hopf equation (|44l) . we obtain the following identities 

S+ -gJ}+^Ra+ =0, (3eC+, 

(45) 

-H--S i?(e-*'^iiior-i«*r)=o, pec-. 

Note that 

/OO 
e-'^'lT/j+^Ra+dt < OO, 
-OO 

and 

e-^P-W^ Rlu.j- ^ o (J^y /3^oo, /3GC-, 

for any n G N. Consequently, we obtain from the leading term asymptotic of both identities (|45l) the 
same result, namely 

1 

MiK,{a)=-— e-'^^'lUj+^Ra+dt. (46) 



The application of the Tauberian type Theorem IA.5I (see Appendix IA.2P leads to 

lim MiK^a) --{-iMiKo), 

a->0+ I 

and, since the matrix AAi is not singular, 

lim KJa) = Kq. 

This proves that the complex stress intensity factor K is continuous in a = 0. 

In order to derive the first order variation of the complex stress intensity factor in the perturbation 
problem, we consider the derivative of l(46|) with respect to a (see Theorem I A. 5p . 

where 



Applying again the Theorem IA.5I we can conclude that 

lim M-iK'Ja) = -iM.2AQ, 

and, finally, the first order variation of the complex stress intensity factor in the perturbation problem 
is given by 



which is consistent with the result l(43|) . 

This asymptotic procedure can be extended to compute the perturbations of high-order terms as 
follows. From the second-term asymptotics of (|45|) we derive the coefficient A* (a), which reads 

■M.A.(a) ^ lim A r te-nur-jmt,,^ ^ ^ _ 



Note that the integral on the right-hand side exists, but the behaviour of the integrand does not allow 
us to take the limit directly. 
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Let us consider the function 



Integration by parts leads to 



where 



9{t,f3)== I -^d^^ (49) 



The integrand in (|49l) is analytic in C and decays exponentially as^-^oo,^eC .As a result 

lim g(t, /3) = 0, /? G C+, 

and therefore the first term in brackets in (l48ll vanishes and we can then write 



27ri 
27ri 
27ri 



/oo 
.g(t,/3)(tIC7]+^i25^+ - Aliii:o)'dt 
-OO 

( [g{t, mm^^R^o - MiK„)] f g'{t, P){tm+^ - MiKo)dt 

L J — oo 

/oo ~iat 

It is now possible to take the limit as /? — > cx) in (|47| . thus we deduce 

M.Ma) ^ lim A e-mur-R^i-M.Ko)^^ 

= -TT- / e-"*(tIC/]+^i?^+ - MiKo)dt. 



We are now in the position to use the Theorem IA.5| so that 

lim Ai2A^(a) = --{-iM2Ao), 
and, since the matrix AI2 is not singular, 

lim Ai,{a) = Aq, 

which proves the continuity of the function A^,{a) in a = 0. 

Moreover, the first order variation of A* (a) can now be derived following the same asymptotic 
procedure presented above and we deduce 

so that 

A'M = (3/2 + ie)Bo, A;'(0) = (3/2 - ie)B*. 
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6 An illustrative example 



In this section, we present an example regarding the computation of the complex stress intensity factor 
K for an interfacial crack loaded by a simple asymmetric force system, as shown in Fig. [31 The loading 
is given by a point force F acting upon the upper crack face at a distance a behind the crack tip and 
two point forces F/2 acting upon the lower crack face at a distance a — b and a + b behind the crack tip. 
In terms of the Dirac delta function 5{-), the loading is then given by 

p+{xi) = -F6{xi+a), 

_ F F 

p {xi) = -—S{xi+a + b)- —6{xi+a-b). 



F/2 



F 


X2 

a 






F/2 


F/2 


1 



2b 



asymmetric loading 



t 


t 




1 i 


F/4 i 


F/4 



symmetric part 



F/2 



F/4 F/4 

anti-symmetric part 



Figure 3: "Three-point" loading of an interfacial crack. 



The loading is self-balanced, in terms of both principal force and moment vectors, and can be divided 
into symmetric and skew-symmetric parts (see Fig. [3]), 

F F F 

{p){xi) = -—6{xi + a) - -^S{xi +a + b) - -^S{xi + a - 6), 

F F 
ipi(xi) = -FS{xi +a) + -^S{xi +a + b) + -^S{xi +a-b). 

For this loading system, the complex stress intensity factor K has been evaluated by means of the 
integral formula l(38|) and the symmetric and skew-symmetric weight functions, HI) and ^ respectively, 
obtaining K = + , where 



(1 + b/a)-^/^-'' + -(1 - b/a)-^^^-' 



K"^ = aF 



cosh(7re)a 



The complex stress intensity factor has been computed for the case = 0.2 and — 0.3 and five 
values of the parameter ry, namely 77 = {—0.99, —0.5, 0, 0.5, 0.99}. The values have been normalized and 
plotted in Fig. [4] as functions of the ratio b/a {a is fixed equal to 1). The symmetric (skew-symmetric) 
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Figure 4: Symmetric and anti-symmetric stress intensity factors as a function of h/a (a is fixed equal to 1) for the case 
v+ = 0.2, V- = 0.3 and: r) = -0.99 (green), Jj = -0.5 (orange), 7] = (red), rj = 0.5 (blue), r) = 0.99 (black). 



stress intensity factor is reported on the left (right) of the figure, the real (imaginary) part is reported 
on the top (bottom). 

Commenting on the results, first, we note from the figure that both K^^ and are not identically 
zero, even though the loading corresponds to a tensile (or opening) mode. This is typical of the interfacial 
crack problem, where there is no symmetry with respect to the plane containing the crack, so that Mode 
I and Mode II are coupled. 

Second, results pertaining to the skew-symmetric stress intensity factor show that Kf and are 
both equal to zero for 6/a = 0, as expected, since in this case the loading is symmetric. As we increase 
b/a and the skew-symmetric part of the loading becomes more and more relevant, the skew-symmetric 
stress intensity factor increases correspondingly. Both the symmetric and skew-symmetric stress intensity 
factors are singular as b/a 1, because a point force is approaching the crack tip. 

The analytical method developed in the present paper allows us also to evaluate the coefficients in 
the higher order terms of asymptotics of physical fields. As an example, we computed the coefficient in 
the second order term. A, which is reported in Fig. [H in analogy with the stress intensity factor in the 
figures above. 

In order to appreciate the magnitude of the skew-symmetric stress intensity factor with respect to 
the magnitude of the symmetric stress intensity factor, the ratio K^/Kf is plotted in Fig. [6] (on the 
left) as a function of b/a. We notice from the figure that the magnitude of may easily reach 20% of 
the magnitude of K^' , a value which is not negligible and has to be taken into account in the view of the 
appHcation of a tensile fracture criterion. Note also that the coefficient A is even more sensitive to the 
anti-symmetric loading, A^ may easily reach the same magnitude as , see Fig. [6] on the right. The 
ratios K^/Kf^ and Aj^/Afj are not reported here since, in this case, they are both equal to the same 
constant, namely K{^/Kf^ = ^n/Afj = a. 
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Figure 5: Symmetric and anti-symmetric parts of A as functions of b/a (a is fixed equal to 1) for the case = 0.2, 
u- = 0.3 and: rj = —0.99 (green), t] = —0.5 (orange), = (red), rj = 0.5 (blue), t] = 0.99 (black). 




Figure 6: Ratios Kf/Kf and Aj*/Af as functions of b/a {a is fixed equal to 1) for the case u+ = 0.2, 
T) = —0.99 (green), r) = —0.5 (orange), 77 = (red), rf = 0.5 (blue), r) = 0.99 (black). 



0.3 and: 



7 Mode III 

The Wiener-Hopf equation relating the Fourier transforms of the singular displacement and the corre- 
sponding traction is given by 

where e — 1^+/^+ + This equation can be immediately factorized as follows 



Assuming that |;73l+ - /3+^^^ and E32 - (3^/^ , as /3 -> 00, /? G C±, we have 

,i/2r— ,4- 6-l-e- 



13. 



1/2 ^32 



0(1), /3^±oo, /3eR, 
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so that, from the Liouville theorem, it follows that both sides of the equation are equal to the same 
constant. 

The weight functions for the Mode III problem are then given by 

_ ^1/2 

[C/3l + (/3)-/3;'/', and ^ 

with corresponding inverse transforms 

1 - i _i/2 f 0, xi > 0, 

-x^ , xi > 0, 



The Betti identity reduces to 

The asymptotics of physical fields are 

Km -1/2 ^III 1/2 3/2s „ + 

(732 = ^=2:1 + —7^x^ + 0(a;/ ), xi 0^, 

V ZTT V ZTT 

^,3j = ^(^±£l^(_,,)l/2 _ 2(5 + eMlII (^^^^3/2 + 0[(_,,)5/2]^ _ 



/27r 3V27r 
with corresponding Fourier transforms 

'^32 = 2 ^+ 4 ^+ '^'^^+ ~* P £ 

H- (l + ^)(& + e)A:„I 3/2 (l-t)(6 + e)A„i 5/2 , ^^^-7/2^ „ _ « ^ ,r- 
IU3JI = ^ /i_ H /i_ + Cl/^- ' ), ~> 00, e L . 

By the same arguments used for Mode I and II, we obtain the formulae 

.0 

Kin = -{l + i) lim / {It/3l(.T'i-xi)(p3)(a:i) + ([/3)(.T'i-a;i)b3l(xi)}dxi, 

A on , ■^ r r i\iTT^( ' ^ d{p3)ixi) , dlp3j{xi) \ 

Aiu = -2{l + t) hm / i lUslix-i^ ~ xi) \- {U3)[Xi - xi) V dxi. 

x[-*o+ I dxi dxi J 

From the solution of the half-plane problem, we get the full representation of weight functions. For 
the lower half-plane: 

^3(/3,.2)^ ^^32^1^1-, 

^3i(/3,a^2) =-^sign(/3)S-2el^l^^ 
a32(/3,X2) =E3"2el''l"^ 

For the upper half-plane the equations are the same subject to replacing |/?| with — and /i_ with 
so that we obtain 
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where 

- AH 
V = : ■ 

/i- + At+ 

The formulae for the evaluation of Km and Am are 
Km 



Am = -2\l- 



- / {(P3>(xi) + |M(2;i)}(-.Ti)-^/^da;i, 

T . w ^ (52) 



dxi 2 dxi J 



The full representation of weight functions after Fourier inversion is, for the lower half-plane, X2 < 0, 

^3(2:1,3:2) = -—^-^ \ {ixi - x^Y^^"^ - i{-ixi - a;2)"^/^| , 

2V7r(/i- + L J 

0'3i(a:^i,a::2) = :r^=7 ; -\i(ixi-X2) ' -{-ixi-X2) ' \, 

cr32(a;i,X2) = - , ^f^^T \ ^ 2:2)"^^^ - i(-ia;i - .T2)"^/n , 

+ 11+) y J 

and, for the upper half-plane, X2 > 0, 

U3ixi,X2) = ^ T {{ixi + X2)~^^^ - ii-ixi + X2y^/^] , 

'J3l{xi,X2) = -—J^^^ (1(1X1 + X2)~^^'^ - {-iXi + X2y^^'^\ , 

|(lXi+X-2)-3/2„^(„j^^^^^2)-3/2|_ 



0-32(3^1, 2:2) - 



We compare now the result obtained for the Mode III with that of Mode I and II. In particular, 
we notice from l(52|) that the skew-symmetric part of Km is proportional to the material parameter 77, 
whereas the skew-symmetric part of the complex stress intensity factor K is proportional to the material 
parameter a, see l(38|) and A comparison of these two parameters is then useful to understand 
whether the skew-symmetric loading is more relevant for Mode III (a > rf) or Mode I and II {r] > a). 
This comparison is shown in Fig. \7\ 




Figure 7: The material parameter « as a function of the material parameter ri for: u+ = 0, u- = 0.5 (green); = 0.2, 
u- = 0.3 (orange); 1/+ = u- = 0.3 (red); 1^+ = 0.5, u- = (blue). 



31 



References 

[I] Antipov, Y.A. (1999) An exact solution of the 3-D-problem of an interface semi-infinite plane crack 
J. Mech. Phys. Solids 47, 1051-1093. 

[2] Bercial-Velez, J. P., Antipov, Y.A. and Movchan, A.B. (2005) High-order asymptotics and perturba- 
tion problems for 3D interfacial cracks. J. Mech. Phys. Solids 53, 1128-1162. 

[3] Bueckner, H.F. (1987) Weight functions and fundamental fields for the penny-shaped and the half- 
plane crack in three-space. Int. J. Solids Struct. 23, 57-93. 

[4] Bueckner, H.F. (1989) Observations on weight functions. Eng. Anal. Bound. Elern. 6, 3-18. 

[5] Erdogan, F. (1963) Stress distribution in non-homogeneous elastic plane with cracks. ASME J. Appl. 
Mech. 30, 232-236. 

[6] Hutchinson, J.W., Mear, M.E. and Rice, J.R. (1987) Crack paralleling an interface between dissimilar 
materials ASME J. Appl. Mech. 54, 828-32. 

[7] Lazarus, V. and Leblond, J.B. (1998) Three-dimensional crack-face weight functions for the semi- 
infinite interface crack - I: Variation of the stress intensity factors due to some small perturbation of 
the crack front J. Mech. Phys. Solids 46, 489-511. 

[8] Mishuris, G.S. and Kuhn, G. (2001) Asymptotic behaviour of the elastic solution near the tip of a 
crack situated at a nonideal interface. Z. Angew. Math. Mech. 81, 811-826. 

[9] Meade, K.P. and Keer, L.M. (1984) On the problem of a pair of point forces appHed to the faces of 
a semi-infinite plane crack. J. Elasticity 14, 3-14. 

[10] Noble, B. (1988) Methods based on the Wiener-Hopf technique for the solution of partial differential 
equations. Chelsea Publishing Company, New York. 

[II] Piccolroaz, A., Mishuris, G. and Movchan, A.B. (2007) Evaluation of the Lazarus-Leblond constants 
in the asymptotic model of the interfacial wavy crack. J. Mech. Phys. Solids 55, 1575-1600. 

[12] Pindra, N., Lazarus, V. and Leblond J.B. (2008) The deformation of the front of a 3D interface 
crack propagating quasistatically in a medium with random fracture properties. J. Mech. Phys. Solids 
56, 1269-1295. 

[13] Rice, J.R. and Sih, G.C. (1965) Plane problems of cracks in dissimilar media ASME J. Appl. Mech. 
32, 418-423. 

[14] WiUiams, M.L. (1959) The stresses around a fault or crack in dissimilar media. Bull. Seism. Soc. 
America 49, 199-204. 

[15] WiUis, J.R. (1971) Fracture mechanics of interfacial cracks. J. Mech. Phys. Solids 19, 353-368. 

[16] WiUis, J.R. and Movchan, A.B. (1995) Dynamic weight functions for a moving crack. L Mode I 
loading. J. Mech. Phys. Solids 43, 319-341. 

[17] Zemanian, A.H. (1968) Generalized Integral Transformations. Interscience PubHshers, New York. 



32 



A APPENDIX 



A.l Material parameters relevant to the interfacial crack problem 

The classical material parameters involved in the problem of the interfacial crack are the bimaterial 
constant, 

J_log.'^+ + (3-4.+ )M- 



27^ ''/n- + (3-4i^-)/z+' 
and the Dundurs parameters 

(see for example Hutchinson, Mear and Rice, 1987). Note that the bimaterial constant e can be written 
in terms of the parameter as follow 

e = — arctanh((i*) — — log — — — p- 

Commonly used are also the parameter 

1 — 1 — z/_ 



b = 



and, mainly for the 3D case, the parameter 

In the study of the asymmetrical case, new bimaterial parameters appear. In particular, the parameter 

^ ~ 2^_(1 - v^) + 2^+(l - v^) ' 
for the plane strain problem (Mode I and II), the parameter 

M- - Ai+ 

for the antiplane shear problem (Mode III), and the parameter 

/== — - — 

for the 3D case. 

All the parameters involved in the problem of the interfacial crack are Hsted in Table EH 
We mention here that another bimaterial parameter, denoted by v, has been introduced by Lazarus 
and Leblond (1998), which is defined by 

(l-J^+)//i+ + (l-J^-)//i- 



(1/m+ + l/Ai-)cosh^(7re) 
and can be expressed in our notations by 
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Table 1: Material parameters relevant to the interfacial crack problem {k± = 3 — 'iu±). 



Symmetric weight function fU} 



Skew-symmetric weight function (U) 



b = 





1 - z/_ 






1 + K+ 


1 + K_ 







ba = 









1 + K_ 






1 - 2i/+ 


1 - 2iy^ 


2/i+ 




— 1 


K_ - 1 

+ — 



1 - 2i/ 1 1 - 2iy^ 



2/i+ 2^_ 

— 1 K_ — 1 



67 



4^4 



4/^4 







AH 




3 — 




4/i+ 



/ = 



3 — K_ 
4a^- 



3 — K I 3 — K_ 



4a*4 



4a*- 



d* = 



6 ^ 2^_(1 - + 2//+(l - 
A*-(^+ - 1) - A*+(^- - 1) 

^-(1 + K+) + /i+(l + K-) 



7* = 



7 ^ A*-(l-2z.+ )+/i+(l-2z._) 
a 2^t_(l - - 2/i+(l - 

- 1) - 1) 

A*-(l + - A«+(l + «-) 



77 = 



= 5 



= ba + f 



As shown by Pindra et al. (2008), this bimaterial parameter has the meaning of "equivalent Poisson 
ratio", in the sense that it satisfies the inequalities < < 1/2, and, when the two materials become 
the same, this quantity reduces to the Poisson ratio of the homogeneous body. 

There are several relations between material parameters that appear frequently in the solution of the 
problem. These relations are listed below for convenience of the reader. 



5/2 



1 + rf* 
1 — (i* 



1/4 
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cosh(7re) = -pr, sinh(7re) = 

cosh ( ^) + sinh (^) = e„, cosh (^) - sinh (^) = 1 
2 1 H~ 



do 1 - d* 



A. 2 Abelian and Tauberian type theorems 



In this section, we will prove a proposition (Theorem I A.sp . which is used in Section 1531 for the purpose 
of evaluation of the perturbation of stress intensity factors and high-order analysis near the tip of the 
interfacial crack. This theorem appears to be a new result, not present in the literature. 

We begin with presenting some known theorems. Although these results can be found in the literature 
(see for example Noble, 1988), they are presented here because of their relevance for the analysed problem. 

Theorem A.l. //, for some a G M, |/(x)| < as x +oo and 

f{x)^x'', x~*Q+, -~l<k<0, (A.l) 

then 



is analytic in = {< e C| Im(t) > a} and 



V'+(t)-r(fc + l)e'^^''+^^t"''~\ t^oo, arg(t)-^ 



<A<^. (A.2) 



Theorem A.2. Let i/j^{t) be a function analytic in Cq, with the property (A . gj) where it is assumed 
that — 1 < fc < 0. Then there exists a function f with support M.^, 



1 



satisfying the asymptotic estimate ^A.l\) . 

Theorem A. 3. Let ijj be a locally integrable function on K. satisfying, for some small e > 0, 

V'(t) = A±|t|-'=-i +0(|tr''"^~'), t^±cx), -l<fc<0. (A.3) 
Then there exists a function f, 

f(x] = 

27r 
satisfying 

fix)^B±\x\'^ + Oi\x\''+'), x^0±, (A.4) 

where 

B± = \ {A+ + sin ^ ^ ^ T ^{A+ - A_) cos ^ ^ ' 



fix) = ;i / V(i)e-"*dt, 



35 



Proof. First we write 

where tpiit) = H{t)'ijj{t) and ^/'2(0 = [1 — H{t)]ij}{t), with H{t) being the unit-step Heaviside function. 
Then, we consider the inverse transform of V'i(i) and, assuming a positive constant a, we deduce 



27r 7o 



J a 



1 

1 /"^ 

=: — / 7/'i(t)e-"*di + 0(l), x-^O, 

where the last equahty is a consequence of the local integrability of the function ipi . From the asymptotics 
l|A.3P i ■ we have ^ A+t~^~^ , t ^ +oo and thus 

h{x) = g{x) + Cl(a;"""{'=+'^'i>), x ^ 0, 

where 



Making the substitution ^ = \x\t, we obtain 



A r°° 

9{x) - / r'-'e— •sn(-)«rf^^ (A.5) 

271" J|a;|a 



and hence 



z:-»o± a; r 27r 



27r 

Considering now the inverse transform of ^^2 {t) , we have 



^{j r'-'cos^d^Ti I r'~'sinCdef> (A.6) 

A+r(-A:) r . 7r(fc + l) 7r(fc + 1) 
sm ^ =F I cos ^ 



f2{x) = — / Mty~"'dt = — / Mty^'dt, 

27r J_oo 27r Jo 



where V'3(0 = i'2{—t) ~ ^-i ^) ^ ^ +oo. From the same argument used above for ,fi{x), we deduce 

/2(x) A_r(-fc) f . ^(fc + i) , . ^(fc + i) ] 

Imi , ,, = < sm ± I cos > . (A. 7) 

It suffices now to add the partial results (|A.6P and (jA.TP to get the final result. □ 
Theorem A. 4. Let Tp be a locally integrable function on M satisfying 

i'it) - A±\t\-\ t^±oo, (A.8) 
and assume that + ^ 0. Then there exists a function f , 



1 r°° 

fix) = — y me-'^'dt, 
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satisfying 



(A.9) 



Proof. The proof proceeds in the same way as for Theorem IA.3| except that the function g{x) in (jA.sp 
is now given by 



1 

-e 



""^^d^ = ^ {- Ci(|x|a) - z sign(a;) [| - Si(|x|a)] } , 



27r 7|2^|a C 



where Ci and Si stand for the cosine and sine integral functions, respectively. From the properties 
Ci ~ log and Si(a;) — *■ 0, as a; ^ 0, we deduce 



A+ 
' 2-K 



\og\xl 



and hence the assertion. □ 

For the purpose of application of these results to the perturbation analysis of the physical fields 
around the interfacial crack, the need arises for an extension of the Theorem lA.2l to the case k This 
extension is provided by the following theorem. 

Theorem A. 5. Let f{x) he the function 

1 /"^ 

I{x) = — y i'+{t)e-"'dt. 

Ifip~^(t) is analytic in C"*" and 

in the closed half-plane — C+ U M, then f{x) — 0, Va; < and 



(A.IO) 



(A.ll) 



lim f(x) = —iai, 

j:-+0+ 



(A.12) 



(A.13) 



Proof. The fact that f{x) =0, Vx < is a direct consequence of the fact that il.'~^{t) is a " + " function. 
Assume now that a; > 0. 

From the assumptions on the behaviour of the function V'^(i)) it follows that ip^(t) ~ ai/t + R{t), 
where tR{t) —* 0, as t ^ oo, t G C (including t ±oo, t e M). 
We write 



/(^) 

The first integral is 



— lim 

2tT a^+oo 



[7/'+(-0e''^* + V+(i)e""*]^^i+ / i'^{t)e-''='dt 



(A.14) 



where 



/ii(a;,a) 



t t 



dt = —2iai 



d^, 
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and 



fi2{x, a)= [i?(-i)e"* + i?(t)e-"*] dt 



-R 



X \x 



so that, taking a = x^^^^, we have fii{x,x^^^^) —inai and /i2(.t, x^^/^) ^ 0, as x 0+. 

Let us denote the second integral in l|A.14p by f2(x,a). Then, using analyticity of in C"*", we 
deduce 

hix,a) = -J^ V'+(Oe-"*dt, 

where Ta ^ {t e C\t = ae'^ <9 <tt}. We write 

f2{x,a) = f2i{x,a) + f22{x,a), 

where 

hi{x,a) = - ^e-"*di and /22(x,a) = -^ R{t) 
Taking again the same a = x^^/^, we obtain 



/2i(x,a;-i/2) = - 



'^1 —jxi 11 

— e dt ~ -ai 



-di — —inai, 



f22{x,a) = - / R{t)e-'-'dt 



as X ^ 0+. 
Finally, 

Jr^„ X \x J 

so that /22(2:,x-i/2) ^ 0, as a: ^ 0+, and we can conclude that 

f{x) —> — (-iTrai - inai) = -iai, 
Ztt 

as a; ^ 0+, which proves l|A.12p . 

Taking the derivative of (|A.10p with respect to x we obtain 

Note that this integral does not exist for real x, so everywhere later we assume Im(.T) > 0. Integrating 
by parts 



fix) 



2^ 



t -I oo 



■'dp 



e-"f dp dt 



2^ 



[til:+ it) ~ ai]' / e-'^Pdp]dt, 



because e "''dp as t ±oo. Integrating again by parts 



[#+W-ai 



dp 



-oo •> -oo 



2^ 



Note that the last formula makes sense also for real x. 



□ 
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An immediate consequence of the Theorem IA.5I is the following 
Corollary. It follows from the Theorem lA.Sl that 

lim f'(x) = — 02. 

Remark. Note that the difference between the Theorems \A.5\ and lA.B is in the assumption about the 
asymptotic behaviour of the transform V'^ as t — > oo. The assumption iA . gj) is weaker, since it is 
equivalent to assume the asymptotic property as lm{t) +00, whereas the assumption IjA . 1 1]) is stronger, 
since the asymptotic property must hold true also as t ^ 00 along the real axis (lm(t) = and Re(t) 
±00). The stronger assumption is needed in order to extend the theorem to the case fc = 0. 

Note also that analyticity of the function 4'^ in the Theorem \A.5\ is essential. To make this clear let 
us consider the following examples. 

In the case k ~ and ^+ + A_ = 0, neither Theorem \A.^ nor Theorem \A.4\ are applicable and the 
only conclusion on the asymptotic behaviour of the inverse transform f{x) is that f{x) is asymptotically 
equal to a constant as x 0, but the value of this constant depends, in general, on the particular 
function ip{t), not only on its asymptotic behaviour as t ±00. For example, if ip{t) = t/{t'^ + 1) then 
f{x) -> iFi/2 as X ^ 0±, whereas if ipit) = (t + l)/{t^ + 1) then f{x) 1/2 =F «/2 as x 0=^. From 
these examples, it is clear that smoothness of the function ip{t) is not helpful. 

Analyticity of the function ■ip{t) = ip^{t) in C^, in addition to conditions f^.3j) guarantees the result 
(Al^) fork = Q and A+ + = 0. 
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